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ABSTRACT Subscripts: 


Analysis of energy emitted from simple or complex 
cavity designs can lead to intricate solutions due to 
nonuniform radiosity and irradiation within a cavity. 

A numerical ray tracing technique was applied to simu- 
late radiation propagating within and from various 
cavity designs. To obtain the energy balance relation- 
ships between Isothermal and nonl sothermal cavity sur- 
faces and space, the computer code NEVADA was utilized 
for its statistical technique applied to numerical ray 
tracing. The analysis method was validated by compar- 
ing results with known theoretical and limiting solu- 
tions, and the electrical resistance network method. 

In general, for non! sothermal cavities the performance 
(apparent emisslvity) is a function of cylinder length- 
to-diameter ratio, surface emisslvity, and cylinder 
surface temperatures. The extent of nonl sothermal con- 
ditions in a cylindrical cavity significantly affects 
the overall cavity performance. Results are presented 
over a wide range of parametric variables for use as a 
possible design reference. 

NOMENCLATURE 

A surface area 

A 0 ut projected surface area at cavity exit 

B radiation interchange factor 

E 0 ut energy emitted per unit time from the cavity 

F geometric configuration factor or blackbody 

view factor 

T absolute temperature 

e surface emisslvity value 

e a apparent emisslvity 


C cylinder 

D disk 

1 , j 1 and j** 1 surfaces 

n total number or nth cavity surfaces 

® space (projected out cavity), infinity 

INTRODUCTION 

The thermal radiation characteristics of partially 
enclosed surfaces (cavities or enclosures) are of great 
interest and have resulted In considerable amounts of 
work in this area. The basic principles of radiation 
exchange for cavity surfaces lead to complex enclosure 
theories. Each cavity surface may emit multiple 
reflecting radiation which can be partially or totally 
absorbed within the cavity or emitted to its surround- 
ings. Depending on the cavity surface properties and 
geometry, the energy leaving the cavity may be composed 
of direct surface emission plus possible reflected 
energy. The term radiosity refers to the rate of 
total radiation (combined emission and reflected) leav- 
ing a surface and the term Irradiation refers to the 
rate of total radiation incident on a surface. The 
nonuniform radiation distribution on cavity surfaces 
results in complex Integral equation formu lations . 
Various theoretical techniques have been applied to 
isothermal and nonl sothermal cavities which generally 
require integration of multiple equarions and numerical 
iteration techniques. The analysis would be further 
complicated by Introducing multiple surface properties 
and various temperature profiles within complex cavity 
configurations. The process of analyzing simple or 
complex cavity geometries may be simplified by applying 
a numerical ray tracing approach to simulate the non- 
uniform radiation propagation. 
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Various types of cavity problems exist. Sparrow 
(1963, 1966) has presented theoretical solutions for 
various isothermal gray-diffuse cylindrical geometries 
with uniform emlssivity throughout the cavity. Due to 
the complicated theoretical solutions required for non- 
isothermal cavities with various internal surface emis- 
slvity values, a ray tracing procedure was applied to 
solve for the energy emitted from the cavities and is 
the subject of this paper. A schematic diagram of the 
analyzed cylindrical cavity is shown in Figs. 1 and 2. 
Figure 1 displays the uniform cylinder wall tempera- 
ture case whereas Fig. 2 displays the linearly varying 
cylinder wall temperature with wall distance case. A 
numerical ray tracing code that simulates the electro- 
magnetic theory of radiation will be used for analyz- 
ing these cavities. The numerical approach solves for 
the total energy out of the cavity from its Individual 
contributing gray-diffuse cavity sections. 


CYLINDER WALL SECTION 
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T D* e D 



CAVITY OPENING 


FIGURE 1. - CIRCULAR CYLINDRICAL CAVITY, CROSS SECTION. 



DISTANCE FROM DISK 

FIGURE 2, - CIRCULAR CYLINDRICAL CAVITY DISPLAYING THE LINEAR WALL 
TEMPERATURE PROFILES. 


There were two basic goals associated with this 
cavity radiation analysis. The first goal was to vali- 
date the numerical ray tracing approach for use in 
studying complex cavities of various geometries, and 
the second was to determine the effects of parametric 
variations on cavity performance. The parametric stud- 
ies would examine the relationships between a range of 
cavity dimensions, surface properties, and temperature 
distributions which may provide insight into radiation 
propagation for nonlsQthermal cylindrical cavities and 
prove useful for design applications. The Isothermal 
analysis of Fig. 1 will be compared to a known theoret- 
ical solution, limiting solutions, and a resistance 
network method as a validation check. 


METHOD OF ANALYSIS 

There are several varieties of thermal analysis 
codes available to model energy flux distributions for 
various geometries. Each code may have different solu- 
tion techniques and advantages. The Net Energy Verifi- 
cation And Determination Analyzer (NEVADA) program 
(Turner, 1988), was selected to simulate the radiation 
propagation within a cylindrical cavity. NEVADA is a 
software package consisting of several programs in 
which a Monte-Carlo mathematical technique is applied 
to radiation propagation. The NEVADA program was 
attractive to use because of its ability to handle dif- 
fuse or specular radiation and, most importantly, com- 
plex geometries. 

In the NEVADA code, a statistical numerical method 
using the Monte-Carlo technique is applied to a ray 
tracing procedure to model radiation exchange. The 
ray tracing procedure mathematically traces emitted 
rays (simulating emitted radiation) as they propagate 
throughout the cavity. Each ray leaving a surface is 
considered a bundle of photons. Each photon bundle 
carries equal, discrete, amounts of energy. The path 
of the bundles (rays) may interact with various sur- 
faces, all of which may reduce the energy of the bun- 
dle. The interacting surfaces may have different 
thermal properties and geometric configurations that 
may affect each bundle's propagation differently. By 
accounting for all the emitted bundles as they propa- 
gate throughout the cavity, percentages of incident 
and absorbed energy at desired locations can be com- 
puted.- The percentages of absorbed energies are then 
applied to energy balance equations. 

ANALYSIS 

In studying cylindrical cavity type configura- 
tions, the overall performance of a cavity is aefined 
as an apparent emlssivity. The apparent emlssivity is 
the actual amount of energy leaving the cavity compared 
with that radiated, by a blackbody at disk temperature 
and projected exit area, or; 



This apparent emlssivity relates direct emitted and 
reflected radiation leaving the cavity to blackbody 
radiation. Due to the nonuniform distribution of radi- 
ation within the cavity (radlosity and i rradi at ion) , 
calculating the energy emitted from the cavity becomes 
complicated. Approximate analytical solutions for the 
radiation exchange Integral equations were first 
derived by Buckley (1927, 1928) and Eckert (1935) ana 
numerically integrated to greater accuracy by Sparrow 
and Albers (1960). 

The numerical ray tracing technique enables one 
to evaluate the radiation heat transfer from a multi- 
sectioned cavity to its various sections and the radia- 
tion projected out the cavity. The numerical ray 
tracing technique actually maps radiation as it propa- 
gates from surface to surface, and can be applied to 
the heat transfer equation as: 

E i->j * e 1 aA 1 B i+j( T t " T j) (2> 



The 1 and j represent various sections of cavity 
surfaces such as Its circular base 0, the cylindrical 
walls C, or Its opening (space). The blackbody view 
factor F Is replaced by the radiation Interchange 
factor B (obtained from the numerl-cal ray tracing 
technique) which represents real surface radiation 
exchange. The radiation Interchange factor Is a func- 
tion of the blackbody view factor and the emlsslvity 
from all energy exchanging surfaces. The radiation 
Interchange factor Is the fraction of emitted energy 
by a real surface 1, that Is absorbed by a real sur- 
face j. Including all reflections from other real sur- 
faces including surface 1. For this cavity analysis 
the surfaces are emitting and reflecting diffuse radia- 
tion. The apparent emlsslvity for a multisection 
cavity where n represents the total number of indi- 
vidual surface sections becomes: 

n 


' T ®) 

1 * | 



Eq. (4) reduces to: 


£ d b d+« 


(9) 


For a limiting solution where gq - cq *) , the cavity 
walls simulate a blackbody absorbing all Incident radi- 
ation making Bq-** - Fq-x® : 


a T- 




( 10 ) 


Thus , for this limiting solution the apparent emis- 
slvlty becomes the disk surface view factor out the 
cavity opening. The uniform cavity wall temperature 
analysis included evaluating the cylindrical walls 
with an emissivity of 0.9 case to simulate near black- 
body effects for a limiting case comparison. In addi- 
tion, various uniform cavity wall temperature cases 
were also solved using the standard resistance network 
method (Holman, 1986), which models only uniform 
radiosity and Irradiation, as shown In Fig. 3. 


By use of the numerical ray tracing approach a 
cavity can be divided into any number of desirable 
sections within computer computational limits. The 
Isothermal and non! sothermal uniform cavity wall tem- 
perature analysis divides the cavity into two sections. 
The circular base 0 and cylindrical walls C com- 
prise the cavity model as shown In Fig. 1. For the 
geometry in Fig. 1, Eq. (3) reduces to: 


e D CTA Q B D-*c' 
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e c aA c B cw 
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(4) 


For the linear cavity wall temperature analysis the 
wall is divided into ten cylindrical sections as shown 
in Fig. 2. 

There are two limiting cases which are used 
in validating the numerical results applied in Eq. (4) 
for the uniform cavity wall temperature analysis. The 
first limiting case applies when the L/R value 
approaches zero. In this case: 
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then B o-*= - F 0~ 


(5) 

( 6 ) 


and Eq. (4) simpl If ies to: 


(7) 


L/R+0 


So for extremely small cylinder cavities where L/R 
approaches zero, the limiting solution for the appar- 
ent emlsslvity becomes the disk surface emissivity 
value . 


The second limiting solution concerns the case 
where the uniform cylindrical wall temperature and the 
environmental temperature approach absolute zero, and 
the surface emlsslvity value of the cylindrical wall 
approaches 1 (complete absorption without emission). 
Thus, if: 


T r ->0 and T +0 
C * 


( 8 ) 


The apparent emissivity results for the isother- 
mal cavities are plotted in Fig. 4 for various L/R 
values and a range of surface emissivity values for 
diffusely reflecting cavity components. For isothermal 
cavities with large surface emissivity values and cyl- 
inder lengths significantly larger than the cavity 
radius, the results simulate an Infinitely deep cavity 
as shown in Fig. 4. As a validation check of the ana- 
lytical method, Fig. 4 displays the apparent emissivity 
distribution for an isothermal numerical ray tracing 
analysis against Sparrow’s theoretical approach and the 
limiting solution of Eq. (7). The numerical ray trac- 
ing results compared exactly with Sparrow’s points and 
agree with the limiting solution for this isothermal 
case. 


The resistance network method was also applied in 
analyzing emitted energy from isothermal cylindrical 
cavity configurations as a possible validation check 
and to evaluate limiting solution criteria for use as 
a possible simplified solution technique. Table 1 dis- 
plays apparent emissivity results from the numerical 
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APPARENT EMITTANCE 



L / R 


FIGURE M. - APPARENT EMITTANCE RESULTS FOR DIFFUSELY RE- 
FLECTING ISOTHERMAL CYLINDRICAL CAVITIES. 


TABLE 1, APPARENT EMISSIVITY RESULTS 
RAY TRACING AND RESISTANCE NETWORK METHOD 
(ISOTHERMAL CYLINDRICAL CAVITY) 



RAY TRACING - NUMERICAL RAY TRACING TECHNIQUE 
R.N.M. ■ RESISTANCE NETWORK METHOD 


ray tracing technique and the resistance network 
method for several cavity dimensions. For relatively 
small cavities (l/R less than 2). for practical pur- 
poses, the resistance network analysis resulted in 
exact agreement with the cavity technique applied with 
Eq, (3) and Sparrow's solution. In these shallow cavi- 
ties the resistance network method accurately predicts 
the energy out the cavity for the entire range of sur- 
face emisslvity values. For cavities with L/R dimen- 
sions larger than 2, the difference between the two 
analysis methods becomes evident. The L/R = 4 
results display a widening discrepancy between results 
at lower cavity emisslvity values. This results from 
the resistance network method requiring assumptions 
for equal energy distributions on individual surface 
sections (radlosity and irradiation). To satisfy the 
requirements the cavity can be sectioned so the energy 
distribution is approximately uniform over each sur- 
face. Once proper cavity sectioning is achieved, the 
results, for practical purposes, agree exactly with 
the numerical ray tracing technique and Sparrow's theo- 
retical results. Therefore, for cavities with L/R 
dimensions greater than 2, the resistance network 
method requires sectioning the cylindrical walls into 
smaller individual surfaces. The problem associated 
with the use of the resistance network method Is prop- 
erly sectioning the cavity to achieve uniform energy 
distribution within each surface, which may also result 
in complicated view factor relationships and a large 
set of simultaneous equations to solve. 

Comparing results derived from the various tech- 
niques applied to mul ti sectioned Isothermal cylindrical 
cavities reveals that using Eq. (3) with the numerical 


ray tracing technique can properly evaluate the total 
energy emitted from the cavities, and the energy emit- 
ted from individual cavity sections. With this vali- 
dation of the ray tracing technique, it can now be 
applied with greater confidence in analyzing more com- 
plicated cavity designs. This type of technique also 
applies to the energy exchanged between Individual sec- 
tions within the cavity. 

APPARENT EMITTANCE 



TEMPERATURE RATIO, CYLINDER TO DISK 
(A) L/R = 0.01. 


APPARENT EMITTANCE 



TEMPERATURE RATIO, CYLINDER TO DISK 
(B) L/R = 0.50. 


APPARENT EMITTANCE 



TEMPERATURE RATIO, CYLINDER TO DISK 
<C) L/R = 1.00. 


FIGURE 5. - APPARENT EMITTANCE AS A FUNCTION OF TEMPERA- 
TURE. 







APPARENT EMITTANCE 



Inspection of the curves reveals that regions exist 
where the apparent emlsslvlty results are dominated by 
cavity cylinder temperatures and surface emissivity 
values, except for considerably smaller cavities. For 
each analyzed cavity dimension, low cylinder wall tem- 
peratures have an insignificant effect on the cavity 
emitted energy, as seen on the left of each Fig. 5 
plot. As the cylinder temperature increases, thereby 
moving right in each Fig. 5 plot, cylinder surface tem- 
peratures and emlsslvlty values have an increasing 
effect on the energy emitted from the cavity. The set 
of figures also include for comparison Sparrow's iso- 
thermal theoretical points and the limiting solution. 
Sparrow's theoretical points are represented at temper- 
ature ratios of one, which corresponds to the Fig. 4 
1 sothermal analysis . 


(D) L/R = 2.00. 


APPARENT EMITTANCE 



(E) L/R - M.00. 
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TEMPERATURE RATIO, CYLINDER TO DISK 
(F) L/R = 8,00. 

FIGURE 5. - CONCLUDED. 

For noni sothermal cavities with uniform cylinder 
wall temperatures, Figs. 5(a) to (f) were developed to 
display the temperature relationships between the disk 
and cylinder wall surfaces, surface emissivity values, 
and L/R effects on the apparent emissivity based on 
Eq. (4), The temperature ratio on the x-axis, cylin- 
der to disk, represents the following equation for 
environmental temperatures at absolute zero: 



The limiting solution in both Fig. 4 and 
Fig. 5(a) at small L/R ratios is the disk surface^ 
emlsslvlty value, Eq. (7). For these extremely shal- 
low noni sothermal cylinder cavities, the cavity wall 
temperatures have little effect on cavity apparent 
emissivity values, as indicated by the nearly horizon- 
tal lines in Fig. 5(a) . 

The limiting solution in Fig. 5(b) to (f) <s the 
disk surface view factor value if the cavity surfaces 
are represented as blackbody surfaces (perfect absor- 
bers of radiation). The view factor values from the 
numerical ray tracing technique were compared to an 
independent solution (Kreith, 1965), which resulted In 
nearly identical values. The limiting solution breaks 
down for nonblackbody surfaces when the L/R dimension 
and the surface emissivity value permit energy (rays) 
to reflect off the cylinder cavity walls and escape 
out the cavity opening before the energy is totally 
absorbed . 

Linear cylinder wall temperature effects (as 
illustrated in Fig. 2) are shown for several cavity 
dimensions with varying surface emissivity values in 
Figs. 6(a) to (d). The major distinguishing factor 
from the previous plotted results is that the tempera- 
ture ratio on the x-axis now represents exit to disk, 
by the following equation for environmental tempera- 
tures at absolute zero: 


Temperature ratio = 



( 12 ) 


APPARENT EMITTANCE 



(A) L/R = 1.00. 


T c - T^ T c 

Temperature ratio - ? t“ 

! D " * 'd 


FIGURE 6. - APPARENT EMITTANCE AS A FUNCTION OF LINEAR 
WALL TEMPERATURES. 
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APPARENT EMITTANCE 



(B) L/R = 2.00. 


APPARENT EMITTANCE 



(C) L/R - *1.00. 

APPARENT EMITTANCE 



(D) L/R = 8.00. 

FIGURE 6. - CONCLUDED. 

Inspection of these curves reveals that the same trends 
occur as In the noni sothermal analysis with uniform 
wall temperatures (Fig. 5), where the apparent emisslv- 
ity is dominated by the cylinder wall temperature pro- 
files. In addition, for large cavity L/R values, the 
plots reveal that upon reaching a specific cylinder 
wall temperature profile, the relationship between sur- 
face emissivity values begins to have a reverse effect 
on cavity performance. The location where the apparent 
emissivity curves converge for the various surface 
emissivity values will be referred to as the cylinder 
emissivity threshold point or just threshold point. 

For practical purposes, the locations of these points 
of convergence are treated as a single point. 


As seen by moving from left to right in the appar- 
ent emissivity plots of Figs. 6(c) and (d), the effect 
of surface emissivity values on the energy emitted out 
the cavity reverses after crossing the temperature 
threshold location. At the left side of the tempera- 
ture threshold location, the low exit surface temcera- 
tures become a minor source of emitted radiation with 
higher surface emissivity values causing less emission 
of disk and cylinder wall energy from the cavity. In 
contrast, to the right of the temperature threshold 
location, the higher exit surface temperatures become 
a major contributing source of emitted radiation where 
a larger surface emissivity value permits an Increase 
of energy emitted from the duct. 

Figs. 7(a) to (d) are also included to display 
and compare cavity results for noni sothermal uniform 
and linear cylinder wall temperatures. Note the tem- 
perature ratio on the x-axis is displayed aoove e or 
the linear and below for the uniform temperature 
curves. Again Sparrow's isothermal theoretical points 
are displayed in the Figures for reference. From this 
comparison, it can be seen that the higher apparent 
emissivity values are ootained for the linear wall tem- 
perature cases than for the uniform wall temperature 
cases. This primarily results from the linear analysis 
containing a higher average cylinder wall temperature 
than the uniform analysis at any given temperature 
ratio, thus producing greater emitted energy. 


TEMPERATURE RATIO, EXIT TO DISK 



(A) L/R = 1.00. 


TEMPERATURE RATIO, EXIT TO DISK 



(B) L/R = 2.00. 


FIGURE 7. - APPARENT EMITTANCE AS A FUNCTION OF UNIFORM 
AND LINEAR WALL TEMPERATURES. 
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TEMPERATURE RATIO. EXIT TO DISK 



(C) L/R = <1.00. 


TEMPERATURE RATIO, EXIT TO DISK 



(D) L/R = 8.00. 
FIGURE 7. - CONCLUDED. 


CONCLUDING REMARKS 

To evaluate simple or complex cavity geometries 
with variable surface properties, a statistical numeri- 
cal ray tracing technique used within the computer code 
NEVADA is an effective tool for obtaining nonuniform 
energy balance relationships between surfaces. This 
method can properly evaluate the total radiant energy 
emitted from complex mu! tl sectioned cavities and the 
radiation exchanged from individual sections within 
the cavity. The cavities need not be sectioned into 
individual surfaces with uniform Incident radiant 
energy distributions, but rather surfaces only to 
define the cavity geometry. For Isothermal cylindri- 
cal cavities the numerical ray tracing technique com- 
pares exactly with other known theoretical solutions 
and limiting solutions. For relatively small cavities 
(L/R less than 2) the resistance network (electrical 
circuit) method accurately predicts the radiation emit- 
ted from the cavity for the entire range of surface 
emissivity values. For larger cavities the resistance 
network requires sectioning the cavity surfaces into 
uniform energy distributed sections which may lead to 
an extremely complicated solution process. 


For noni sothermal cylindrical cavities with either 
uniform or linearly varying cylinder wall temperatures , 
the energy emitted from the cavity is a function of 
cavity surface emissivity value and both cavity size 
and surface temperature except for relatively snaTcw 
cavities. For such cavities, with cylinder wall 'ength 
smaller than cavity radius, emitted energy becomes 
much less dependent on cylinder surface temperatures 
and more dependent on disk surface emissivity values. 

For cavities with linear cylindrical wall temperature 
profiles, there exist temperature threshold locations 
which cause a reversal of the surface emissivity effect 
on the energy emitted out the cavity. 

By plotting the theoretical points with the numer- 
ical ray tracing results and with limiting solutions, 
one can accurately grasp the effect of surface tempera- 
tures and surface emissivity values on energy emitted 
from various cavity configurations. This numerical 
ray tracing technique can enable one to analyze the 
nonuniform energy distributions within cavities which 
are separated Into various sections with different Tem- 
perature profiles and different emissivity values. 
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